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With the improvement in fabrication techniques it is now possible to produce atom-like semi- 
conductor structures with unique electronic properties. This makes possible periodic arrays of 
nano-structures in which the Coulomb interaction, polarizability, and tunneling may all be varied. 
We study the collective properties of 2D arrays and 3D face centered cubic lattices of singly-charged 
nano-spherical shells, sometimes called 'quantum-dot quantum wells' or 'core-shell quantum dots.' 
We find that for square arrays, the classical groundstate is an Ising anti-ferroelectret (AFE), while 
the quantum groundstate undergoes a transition from a uniform state to an AFE. The triangular 
lattice, in contrast, displays properties characteristic of frustration. Three dimensional face-centered 
cubic lattices polarize in planes, with each layer alternating in direction. We discuss the possible 
experimental signals of these transitions. 
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I. INTRODUCTION 

With improvements in the colloidal fabrication of 
structures on the micro- and nanometer scale, new ar- 
eas of study have been made possible. In particular, re- 
searchers are now able to produce atom-like electronic 
devices with their own unique spectra and shell struc- 
tures P, Q . The distinct characteristics of such devices 
may be controlled through fabrication, in turn producing 
structures with unique electronic properties. 

Periodic arrays of atoms have been studied exten- 
sively for nearly a century. However, we may now ex- 
amine periodic arrays of nano-structures in which the 
tunneling, polarizability, and Coulomb interaction may 
all be varied, thus giving us the power to control an elec- 
tronic wavefunction to a degree without an atomic ana- 
logue. Small multi-dot systems show coherent phenom- 
ena in experiments 3] while disordered arrays of colloidal 
quantum dots display variable long-range hopping. ^5] It 
should soon be possible to examine collective electronic 
excitations of periodic nanostructures. With this in 
mind, we choose a nano-structure which will confine the 
electron on the nanoscale in such a way as to bring about 
novel characteristics. 

In this paper we consider electrons confined to spher- 
ical shells, sometimes called "quantum dot quantum 
wells" (QDQW) or "core-shell quantum dots." QDQWs 
are heterostructures where layers of different semi- 
conducting materials alternate in a single nano-crystal. 
One well characterized example is the CdS/CdSe/CdS 
system, i.e. a CdS core and outer shell with a CdSe in- 
ner shell Q. In this case each layer is on the order of 
1-2 nm and is separated by a large electronic band gap, 
thus allowing the investigation of quantum confinement 
in a geometry in which electrons occupy the surface of a 
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sphere. CdS/CdSe/CdS QDQW's have also been found 
to have an electron-g factor which varies as a function of 
quantum well width as well as a transverse spin lifetime of 
several nano-seconds at temperatures approaching room 
temperature [8]. 

We choose QDQW's for two reasons. First, be- 
cause they have already be fabricated and are well- 
characterized, and second, because their excluded core 
structure yields properties not found in regular quan- 
tum dot systems. Previous research has found that 
in a ID array of singly-charged nano-rings, a quantum 
phase transition (T — 0) occurs which allows the system 
to spontaneously break symmetry into anti-ferroeletric 
(AFE) alignment. Studies of the 2D ring problem re- 
vealed another phase transition from random orientation 
to a striped AFE state, analogous to unitary transition 
of the 2D XY model. The focus of this paper is to show 
an analogous transition occurs in an ordered array of 
QDQW's. 

In the section below we describe our model, and then 
our methods for analyzing it in the classical and quantum 
mechanical limits. We then present results for the 2D 
square lattice, the 2D triangular lattice and the 3D face- 
centered cubic (FCC) lattice. Both the 2D square lattice 
and 3D face-centered cubic lattice systems yield a phase 
transition from random (classical) or uniform (quantum 
mechanical) distribution to an anti-ferroelectric state. 



II. THE MODEL 

In our model, each QDQW is treated as a infinites- 
imally thick spherical shell of radius R with an electric 
charge on the surface. Charging might be achieved either 
by doping, tunneling from a backgateor be adjusting an 
electrochemical potential. The radial degree of freedom 
can be neglected when the gap between radial excitations 
is large compared to that of orbital eigenstates. We then 
consider periodic arrays of singly-charged QDQW's in ei- 
ther two or three dimensions and define a distance, D, 
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between neighbors. In 2D we choose D as the lattice con- 
stant so that D > 2R. In a 3D face-centered cubic lattice 
D is half the distance between vertices of the FCC cells 
so that the separation between nearest neighbors is given 
by \f2D. We then consider the electrostatic interaction 
solely between nearest neighbors. Though the Coulomb 
repulsion acts over large distances, we assume that there 
is sufficient screening so that only nearest neighbor inter- 
actions are significant. Furthermore, the insulating shell 
around each QDQW allows us to ignore electron tunnel- 
ing between dots. 



A. Classical Analysis 

Before solving a quantum mechanical problem it is of- 
ten helpful to look at the classical case, which is usually 
easier to solve and provides an guide to the quantum 
ground state. Classically the electron is a point charge 
constrained to the surface of a sphere. Because of this 
constraint, it has a dipole moment (with respect to the 
shell center) of constant magnitude but random direc- 
tion. Adjacent moments interact leading to the polariza- 
tion instability, as we will show below. 

There is only one energy scale in the classical problem, 
the Coulomb energy i?c = /D. The electrostatic energy 
of the array is given by 
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FIG. 1: A schematic picture of the groundstate of classical 
point electrons in a 2D array of spherical shells. Each shell 
is of radius R and separated from its neighbor by a distance 
D. The 2D ordering is antiferroelectric for the infinite system 
size. 

B. Quantum Mechanical Analysis 

Instead of a single point charge on a sphere, in the 
quantum mechanical problem the electron is treated as a 
charge distribution constrained to the surface of a spher- 
ical shell. We now have two competing energy terms: 
the kinetic, which aims to spread the wave function, and 
the potential, which is strictly the Coulomb interaction 
between nearest neighbors. Our starting Hamiltonian is 
H = T + U where the kinetic energy is given by: 
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where is the location of the i-th electron, i and j are 
restricted to nearest neighbors, and the factor of 1/2 cor- 
rects for double counting. We define the expansion pa- 
rameter e = R/D and expand the potential to second 
order: 
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where Ec = / D is a measure of the Coulomb interac- 
tion between shells. Here we specify the position of each 
charge by a unit vector Si pointing from the center of 
the z-th sphere to the charge on that sphere. The unit 
vector Di j points from the center of the i-th to the cen- 
ter of the j-th sphere. We have also assumed that the 
lattice possesses reflection symmetry. From eq. ([2]) we see 
that the system has an AFE Heisenberg interaction and a 
symmetry breaking term. The symmetry breaking term 
encourages alignment of adjacent spins in the direction 
of the lattice vector connecting them. Further analysis 
requires knowledge of the lattice vectors of the specific 
system. 



where Eq = h^/2m*R^, and m* is the effective mass of 
the electron. In moving to angular variables, we write 
the unit vectors in eq® in terms of angle variables, Si — 
sin 9i cos X + sin 9i sin (f>iy + cos 9i z. If we measure our 
energy in units of Eg, our scaled Schrodinger equation is 
then 
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where A is the ratio between the Coloumb and confine- 
ment energies, A — E^/Eg and E is the dimensionless 
energy. 

We calculate the groundstate using a variational ap- 
proximation. For the z-th ring we choose a variational 
wave function, using only the first four spherical har- 
monics (Yim): 

'0var(6'j,0j;Q;j;/3j,7i) = cos(/3i)ybo -I- sin(;3i) cos(ai)yio 



-t- sin(/3^) sin(ai)y ^ sin(6'i) cos((/)i - 7.,)(6) 

where we have used time reversal symmetry to force the 
wave function to be real. We are able to use this form 
since at low energies, each electron will occupy one of the 
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FIG. 2: Spherical plot showing the charge polarization of a 
single lattice site for the 2D square array. Left is a uniform 
surface charge density (I/Att) represented by a mesh sphere 
of radius 1 + l/47r surrounding a solid sphere of unit radius. 
Right is a plot of the polarized charge density, again repre- 
sented by a spherical mesh with a radius 1 + surrounding 
a solid sphere of unit radius. The charge is not strongly lo- 
calized since we are dealing with only the lowest spherical 
harmonics. Nearest neighbors display the same behavior but 
with opposite polarization. 
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FIG. 3: A plot of the classical phase transition from random 
alignment to the anti-ferroelectric (APE) state Fig. ([II A|) 
through 7 = l3Ec, a measure of interaction strength, at zero 
temperature. Here 5 is a measure of the system's AFE order. 
Systems are square with periodic boundaries of size N=10, 
N=20, and N=30 and have R/D = 0.417. Larger system 
sizes have sharper transitions. 



lowest energy states. Even at full polarization the charge Cv 
is not strongly localized; this is demonstrated in fig.© If 12 

the system will not polarize in this approximation, it will 
not be able to do so if we include higher energy states. 
This is done for a unit cell that respects the periodicity 0.08 
found in the classical simulations. Separation between 
neighbors is assumed to be sufficiently large to neglect 
tunneling effects and enable the use of the Hartree ap- 0.04 ^ n = 20 

proximation. ^ 
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III. 2D SQUARE LATTICE 2 4 6 8 10 



We first analyze the 2D square lattice assuming peri- 
odic boundary conditions. Below we present the classical 
and quantum mechanical analyses. 

A. Classical Phase Transition 



FIG. 4: A plot of the heat capacity when R/D = 0.417 for 
N=10, N=20, and N=30 against 7 = fiEc, a measure of inter- 
action strength. Here Cv is the system's heat capacity. For 
larger system sizes, the peak around the transition point is 
sharper. 



The nature of the symmetry breaking term in the 
square lattice can be clarified by resolving the interac- 
tion into components: 
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The first term is an irrelevant constant; the second 
term breaks the rotational symmetry, favoring spins 
aligned along the z-axis. The third term is a Heisenberg 
/erroelectric coupling while the last and larger term is 
an Ising anMerroelectric coupling. An AFE state, when 



possible, will maximize the second and fourth terms, min- 
imizing the total energy. This is confirmed by Monte 
Carlo analysis. 

Monte Carlo 6j simulations were performed on 10x10, 
20x20 and 30x30 arrays. The low temperature 
[e^/kTD) » 1) configuration is an antiferroelectric 
(AFE) pattern. Fig. ([Ij. Furthermore, a phase tran- 
sition from the unordered state to an AFE aligned state 
was found numerically. The transition is Ising-like, with 
a discontinuous change in the staggered polarization and 
a peak in the specific heat which is expected to sharpen 
into a divergence as the system size is increased. The 
value of the temperature at the peak extrapolates to 
kTc — XXX for an infinite size system, which should 
be compared to the expected value of 2.24 for the 2D 
Ising model. |10l| 
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FIG. 5: A plot of two of the eigenfrequencies of the AFE 
square lattice as a function of the in plane wave- vector, k' . 
The upper surface has been displaced upward by one unit 
for clarity. Without the displacement the two surfaces touch. 
The other two normal modes are equivalent to those plotted, 
but rotated by tt/2. Note that all modes are gapped. 



We can expand about the minimum energy state in 
order to calculate the normal modes. Since the polariza- 
tion has a checkerboard pattern, our basis is rotated by 
7r/4 vifith respect to the original axis, and we introduce 
x' — (x + y)/ and y' — {x — y)/ \f2. Our unit cell has 
two elements, an "up" site and a "down" site. Denote 
the position of a point charge on an "up" site by 



= uf^ X + U^^^H + ^l-(«f))2 + (u(^))2 i (8) 

while that on a "down" site is: 

= ) i -f d(^) y - yi-(df))2 + (d(^))2z (9) 
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We then solve coupled differential equations 
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working to second order in the displacements. We find 
find four normal modes; the first two have frequencies: 
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and correspond to in-phase/out-of- phase displacements 
in the x' direction. These frequencies are plotted in figure 
([5]) as a function of k' . The second two normal modes 
(not plotted) are equivalent to the first but rotated by 
tt/2, with x' y' . Note that all are gapped. 



B. Quantum Phase Transition 

It is not a priori obvious what will happen in the 
quantum mechanical case at zero temperature, since the 



quantum confinement introduces a new energy scale. 
Rather than solve for the variational groundstate of an 
N X N system directly, it was assumed that the quan- 
tum groundstate would possess the same symmetry as 
the classical system, so that only a 2 x 2 unit cell of an 
(assumed) infinite system was analyzed. This has the 
advantage of following the correspondence principle as 
A — > oo . A plot of the variational energy for e = ^ , re- 
veals that at small values of the coupling constant. A, 
a single minimum exists. Fig. ©. This minimum cor- 
responds to /3 — 0, i.e. a uniform distribution for the 
electronic wavefunction. As A increases past a critical 
value, two new minima emerge at /? = ±-| indicating 
a quantum phase transition. These values correspond 
to charge localizing normal to the plane. Furthermore, 
nearest neighbors polarize in opposite directions. 

With the simple form of i/'var of eq. (jG)) and the approxi- 
mate interaction of eqlZl we can solve analytically for the 
minima of the variational energy and find: 
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where a = tt/2 and A = 0, and the critical value for the 
interaction is given by Ac = (22i44-3i57r^)e^ ~ 1.756/e^. 
Plotting (3 for several values of e reveals the phase tran- 
sition. Fig. ([7]). Here, the analytical minimum result is 
confirmed numerically. 

As e shrinks, the transition becomes weaker, requiring 
the coupling constant, A, to be large for the transition to 
occur. This phenomenon emerges from the geometry of 
the array. When the separation between shells becomes 
larger, the effective transition radius of each shell also 
increases. Fig. ([5]). The maximum value of e can have is 
1/2, when the spheres touch. When e is at its maximum 
there is a minimum R below which quantum mechanics 
dominates and the system does not polarize: i?min ~ 
7aom/m* . Thus we have the curious result that very 
small shells do not spontaneously polarize. 



IV. 2D TRIANGULAR LATTICE 

Following our analysis of the 2D square lattice, we 
ues to the 2D triangular lattice. Experimental arrays 
made from the self-assembly of colloidal particles would 
likely stack in triangular arrays, rather than square 
ones. For the ring system the triangular lattice sponta- 
neously breaks symmetry into a striped AFE phase with 
a six-fold degenerate groundstate. However, in that case 
the spins order in the plane, not perpendicular to it. 

For a pure AFE Ising system the triangular lattice is 
frustrated. JJLj While eq.© is not a simple Ising AFE 
interaction, it is clear that frustration will also be present. 
Monte Carlo simulation of eqn.([2] ) found multiple low 
energy states indicative of a glassy phase. There was no 
peak in the specific heat nor any abrupt change in the 
polarization to indicate an ordered phase. This meant 
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FIG. 6: A plot of the energy functional from the quantum 
variational calculation for the 2D square lattice. Here, e = 5, 
thus the ratio of shell radius to separation is at a maximum 
(the close-packed condition). A single minimum for small 
values of the coupling constant, A, occurs at /3 = 0, where f5 is 
the mixing parameter found in the variational wavefunction. 
This state corresponds with a uniform charge distribution. 
As A increases, two minima emerge at /3 = These values 
correspond to charge localizing normal to the plane of the 
lattice. 




FIG. 7: A plot of the quantum phase transition in the mix- 
ing parameter, /3, of the variational wavefunction through the 
coupling constant, A, for several values of e. This represents 
the system's transition from uniform distribution to antifer- 
roelectric order, through charge localization. Fig. ((2}. Ana- 
lytical values of /? are given by Eq. 1131 Here we see a good 
fit between numerical and analytical solutions. 



there was no basis for a quantum mechanical ansatz to 
the ground state. 



V. 3D FACE-CENTERED CUBIC LATTICE 

The next logical progression is to stack 2D triangu- 
lar arrays forming a face-centered cubic (FCC) lattice. 



FIG. 8: The critical value of the coupling constant A, plot- 
ted against the ratio of shell radius to shell separation, e. 
It is at these values of A that the transition in the system 
occurs, given e. As separation increases, the coulomb interac- 
tion grows weaker, thus an even weaker confinement energy 
is required for the transition to occur. This corresponds with 
a larger physical radius for each shell, A S> 1. 




FIG. 9: A typical low temperature state that anti-ferroelectric 
(AFE) in alternating horizontal planes. Each triangular lat- 
tice row has a common polarization in the plane, while neigh- 
boring rows have opposite polarization in the plane. 

We chose this model since, in the close-packed limit, we 
expect it to be more physically reproducible through fab- 
rication. Each lattice point has twelve nearest neighbors: 
four above, four in the plane, and four below. Again, the 
non-nearest neighbor interaction is considered negligible. 



A. Classical Analysis 

Through numerical Monte Carlo simulation, a mini- 
mum energy was shown to exist when the system is in 
an AFE state by row. That is, each triangular lattice 
row of the FCC structure has opposite polarization. Fig. 
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FIG. 10: A plot of the classical phase transition from ran- 
dom alignment to the anti-ferroelectric (AFE) state Fig. (|V]) 
through 7 = l3Ec, a measure of interaction strength, at zero 
temperature. Here 5 is a measure of the system's AFE order. 
Systems are face-centered cubic with periodic boundaries of 
size N=4, N=6, and N=8. Larger system sizes have sharper 
transitions. 



(fV)) . Furthermore, a phase transition occurs from the 
unordered state to the AFE aligned state. 

The temperature of the transition is determined by 
the interaction strength found in the Boltzmann factor, 
7 = f3Ec — jSe^ /D, where e is the charge of an electron, 
\/2D is the shell separation, and /3 = 1/kT. One can see 
that with increased system size comes a sharper transi- 
tion to the AFE state, Fig. (jV A[) . Numerical simulations 
indicated a persistent wobble about the in-plane polar- 
ization. Analytic expansion of the interaction about the 
uniform unpolarized state indicated an instability that 
pushes the polarization slightly off-axis. An expansion 
about the new minimum indicated restoring forces that 
were only quartic in the deviation from this state. 



B. Quantum Mechanical Analysis 

To confirm our results from classical analysis, we turn 
to the Schrodinger variational principle as described ear- 
lier. We assumed that the low energy state in our quan- 
tum mechanical computations was polarized in the plane 
and alternated from one layer to the next, ignoring the 
slight wobble found in the classical analysis. (This will 
only give us an upper bound on the phase transition 
boundary). After completing the variational calculation, 
we once again find the charge to localize on each sphere 
to give the array an anti-ferroelectric alignment by tri- 
angular lattice row. Fig. (jV])- Furthermore, a quantum 
phase transition is found from uniform alignment to this 
state through the mixing parameter, /?, found in the vari- 
ational wavefunction. Fig. (fT2)) . 

A plot of the energy functional under the close-packed 
condition, e = reveals that at small values of the 
coupling constant. A, a single minimum exists. Fig. (1111) . 




FIG. 11: A plot of the energy functional from the quan- 
tum variational calculation for the 2D square lattice. Here, 
e — thus the ratio of shell radius to separation is at a 
maximum (the close-packed condition). A single minimum 
for small values of the coupling constant, A, occurs at /3 = 0, 
where /3 is the mixing parameter found in the variational 
wavefunction. This state corresponds with a uniform charge 
distribution. As A increases, two minima emerge at P — 
These values correspond to charge localizing along the plane 
of the lattice. 
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FIG. 12: A plot of the quantum phase transition in the mix- 
ing parameter, /3, of the variational wavefunction through the 
coupling constant, A for several values of e. This represents 
the system's transition from uniform distribution to antiferro- 
electric order, through charge localization. Analytical values 
of P are given by Eq. 1141 Here we see a good fit between 
numerical and analytical solutions. 



This minimum corresponds to /3 = 0, i.e. a uniform 
distribution for the electronic wavefunction. As A in- 
creases past a critical value, two new minima emerge at 
(3 — ±7r/4. These values correspond with charge localiz- 
ing along one of the axis of the shell. Furthermore, each 
triangular lattice row polarizes along opposite axes. 
Analytically solving the energy functional for its min- 
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FIG. 13: The critical value of the coupling constant A, plot- 
ted against the ratio of shell radius to shell separation, e. 
It is at these values of A that the transition in the system 
occurs, given e. As separation increases, the coulomb interac- 
tion grows weaker, thus an even weaker confinement energy 
is required for the transition to occur. This corresponds with 
a larger physical radius for each shell, A ^ 1. 

ima, we find the ground state values of the variational 
parameters: 

riarccos(^) if A > A,; 

\0 if A < Ac. ^ ' 

and where a = tt/2 and A = 7r/2, where the critical value 
for the interaction is given by Ac = ^ ■ Plotting /3 for 
several values of e reveals the aforementioned phase tran- 
sition. Fig. (|12|) . Here again, the analytical minimization 
result is confirmed numerically. As e shrinks, the transi- 
tion becomes weaker, requiring the coupling constant. A, 
to be large for the transition to occur. This phenomenon 
emerges from the geometry of the array. When the sep- 
aration between shells becomes larger, the effective tran- 
sition radius of each shell also increases. Fig. 

It should be noted that this transition is weaker than 



that discovered for the 2D square lattice. Not only is 
the minimum coupling constant larger for each transition 
to occur, but also Ac increases at a quicker rate with 
increased separation. Fig. ([T^. 



VI. CONCLUSIONS 

For a 2D square array and 3D FCC lattice of singly- 
charge quantum-dot quantum wells, a phase transition 
from the unordered state to a anti-ferroelectric state was 
found to occur at low energies through classical Monte 
Carlo simulation. Furthermore, this transition is con- 
firmed through a quantum mechanical treatment of the 
system in which only nearest neighbor interactions are 
considered. The transition is a product of the geometry 
of the system, specifically dependent upon the ratio of 
shell radius to shell separation. Finally, this transition is 
said to be "Ising-like," which opens the doors to a host 
of known collective system behaviors. 

For the square lattice we know that in order to 
have polarization, A > 1.75/e^, which means that 
R^m* /aoD^ ,m > 0.875, where qq is the Bohr radius. 
For CdSe the effecive mass is m*/m = 0.13, a typical 
layer thickness is 0.43nm, and the capping (insulating) 
layer i about l.Gnm.Q If we assume that the shell is two 
monolayers thick, then we need shells with a radius of 
about 7nm in order to have a polarized ground state. If 
the shells have a smaller radius the kinetic cost to localize 
the electrons is too high. Charging the shells may prove 
to be too difficult at first; a simpler experiment would be 
to optically excite electron-hole pairs in a transverse elec- 
tric field so as to create "chains" of elecron-holes across 
the core-shell array. 
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